Linear plasmon dispersion in single-wall carbon nanotubes 
and the collective excitation spectrum of graphene 
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We have measured a strictly linear n plasmon dispersion along the axis of individualized single wall carbon 
nanotubes, which is completely different from plasmon dispersions of graphite or bundled single wall carbon 
nanotubes. Comparative ab initio studies on graphene based systems allow us to reproduce the different disper- 
sions. This suggests that individualized nanotubes provide viable experimental access to collective electronic 
excitations of graphene, and it validates the use of graphene to understand electronic excitations of carbon nan- 
otubes. In particular, the calculations reveal that local field effects (LFE) cause a mixing of electronic transitions, 
including the 'Dirac cone', resulting in the observed linear dispersion. 

PACS numbers: 73.20.Mf,73.22.-f,78.20.Bh 



Single-wall carbon nanotubes (SWNT) and its parent com- 
pound graphene are archetypes of low dimensional sys- 
tems with strongly anisotropic and unique electronic prop- 
erties which make them interesting for both fundamental re- 
search and as building blocks in nanoelectronic applications 
Their electronic bandstructure is frequently studied. In 
graphene, the linear band dispersion at the Fermi level, the 
'Dirac cone', leads to unique characteristics in nanoelec- 
tronic devices 13] ■ One can expect a strong analogy between 
graphene and isolated SWNT for excitations along the sheet 
and along the tube axis, respectively. Within the zone-folding 
model, i.e. neglecting curvature effects, the graphene band- 
structure is sliced along parallel lines when the sheet is rolled 
up into a cylinder. The result are characteristic van Hove sin- 
gularities (VHS) in the density of states (DOS)|[l. Bulk (i.e. 
bundled) SWNT show an optical absorption peak at ~ 4.5 cV 
due to transitions of the tt electrons In vertically aligned 
SWNT (VA-SWNT) one finds the same peak position for on- 
axis polarization and an additional peak for perpendicular po- 
larization at ~ 5.2 eV lH. Further information can be ob- 
tained from collective electronic excitations (plasmons) be- 
yond the optical limit |0] (i.e. momentum transfer q > 0). 
Angle resolved electron energy loss spectroscopy (EELS) as- 
sesses the detailed plasmon dispersion QHtl, but it is so far 
missing for freestanding isolated sp^ carbon systems. 

Models based on the homogeneous electron gas @], or the 
tight-binding scheme ifloifTTIl have been used to describe these 
excitations. The former are however bound to metallic sys- 
tems. The latter have provided valuable insight and predic- 
tions for the properties of isolated sheets, tubes, and assem- 
blies of these objects; in particular, they have predicted an 



almost linear plasmon dispersion for isolated systems. How- 
ever, the tight binding results neglect screening beyond the tt 
bands, and they depend on parameters that hide the underly- 
ing complexity. No realistic parameter-free calculations have 
been performed to predict the plasmon dispersion in these 
systems, nor has its origin been analyzed. Instead, ab ini- 
tio spectroscopy calculations have dealt with absorption spec- 
tra (q 0) for SWNT E [H H, and plasmon disper- 
sions in graphite ifisl [Till : most other available calculations 
are ground-state or bandstructure ones ITtIi . The prediction, 
comparison and interpretation of the full dispersive electronic 
excitations of isolated SWNT and graphene sheets calls for 
new experiments and for ab initio theoretical support going 
beyond bandstructure calculations. 

Indeed, electronic excitations imply a self-consistent re- 
sponse of the entire system and have therefore to be described 
in terms of bandstructure and induced potentials. In solids, the 
latter consist of microscopic induced components (local field 
effects l Usll ) and a macroscopic induced component. The lat- 
ter is responsible for the difference between interband transi- 
tions as measured in absorption, and plasmons as measured in 
loss spectroscopies lfl9ll . In isolated systems, the macroscopic 
component will naturally vanish. Absorption peaks and loss 
peaks atq^O will hence coincide. The microscopic part, in- 
stead, can become dramatically important when macroscopic 
screening is low. Therefore, theory based on bandstructure 
alone will be unable to describe isolated systems, whereas ex- 
periments performed on bundles or graphite are not represen- 
tative for isolated SWNT and graphene. There is hence an 
important gap in our understanding of the properties of these 
isolated sp^ carbon systems. 
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The present work is meant to bridge this gap with a detailed 
EELS study on freestanding mats of VA-SWNT, and corre- 
sponding ab initio calculations for graphene sheets. Our stud- 
ies allow us to give answers to several important questions, 
by ( i) distinguishing a localized perpendicular and a strictly 
linear on-axis plasmon dispersion in isolated single-wall nan- 
otubes; (ii) showing the quantitative similarity to electronic 
excitations in graphene; ( Hi) analyzing the impact of local 
field effects on the linear plasmon dispersion, in particular the 
mixing-in of low-energy transitions; (iv) quantifying the im- 
portance of interactions between neighboring sheets or tubes. 

2 |xm to 7 |J.m thick VA-SWNT material was directly 
grown by catalytic decomposition of alcohol and subsequently 
detached from the supporting silicon wavers by floating off in 
hot water and transferred onto Cu grids lEoll . The nematic 
order as well as optical properties 



2 in and local morphol- 



ogy in transmission electron microscopy of VA-SWNT Il22ll 
have been studied earlier. The VA-SWNT are aligned within 
25° and typically packed in small bundles with less than 10 
nanotubes, each with a diameter of about 2 nm. The angle 
resolved loss function of the VA-SWNT was measured in a 
purpose built EELS spectrometer ll23ll . Earlier comparative 
EELS studies were performed on a cleaved single crystal llSIl 
of graphite or bundled and magnetically aligned SWNT Jsl]. 
In the present study we set an energy and momentum resolu- 
tion of 200 mcV and 0.05 A"^ at a primary incidence energy 
of 172 kcV. Concerning our ab initio simulations for isolated 
single- and double-layer graphene, we start from DFT-LDA 
ground state calculations [24], using a plane-wave basis set ' 
and norm-conserving pseudopotentials of Troullier-Martins 
type 12511 . The loss function is determined within the random 
phase approximation (RPA) using the DP-code The local 
field effects (LEE) that originate from the induced Hartree 
potential are taken into account, comprising in-plane local 
fields with spatial variations on the atomic scale (Umklapp 
effects). Moreover, all contributing valence-electron bands (tt 
and a as well as empty states) are included. Moving on to 
our EELS measurements on VA-SWNT, we first inspect the 
loss function at the lowest momentum transfer (0.1 A"^), de- 
picted topmost in Fig. la. We observe peaks corresponding to 
the TT and the more structured tt + cr plasmon at 5.1 eV and 
17.6 eV, respectively. These values are remarkably low for 
sp^ hybridized carbon which indicates minimal macro- 

scopic screening and therefore fits well to the morphological 
investigations of Einarsson et al. 12211 . TEM micrographs in 
Fig. lb(&c show the cross section and side view of this thin 
bundled VA-SWNT. For increasing momentum transfer, the 
loss spectra change significantly (shown up to 1.0 A"^). Ap- 
parently, both, the TT and tt + cr plasmon split into two distinct 
contributions; one is localized like in a molecule and another 
is dispersive like in a solid. We interpret the localized re- 




' We used 3481 k-points and an energy cutoff of 28 Hartree. 
^ http://www.dp-code.org; V. Olevano, et al, unpublished. 
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FIG. 1: (a) measured loss function of freestanding VA-SWNT at 
equidistant q from 0.1 A"^ (top) to 1.0 A"^ (bottom). And TEM 
micrographs of the cross section (b) and side view (c) of the thin 
bundled VA-SWNT 



sponse as a contribution with momenta perpendicular to the 
axis, where an isolated quantum wire becomes a confined dot. 
The dispersive response belongs to momenta parallel to the 
nanotube axis, where we encounter an isolated one dimen- 
sional wire. The simultaneous, but distinct, observation of 
these two intrinsic aspects of nanotubes - being solid-like at 
the one hand and molecule-like at the other - is a fingerprint 
of fully individualized viz. ideal wires in terms of plasmon 
excitations. Following this interpretation, the extrapolation 
of the corresponding tt plasmon position to the optical limit 
{q — > 0) predicts values of 4.6 cV and 5.1 eV for the on-axis 
and perpendicular component, respectively (see Fig. 2). This 
is in excellent agreement with the optical absorption findings 
of Murakami et al. jst]. Focusing on the dispersive on-axis 
part of the tt plasmon we find a strikingly linear behavior 
(filled diamonds in Fig. 2d) up to about 1 A"\ covering more 
than one third of the Brillouin zone. The linear behavior is 
totally different from the parabolic tt plasmon dispersion of 
bulk graphite (filled triangles in Fig. 2b), while bulk aligned 
SWNT (filled circles ^ in Fig. 2c) represent an intermedi- 
ate situation. Therefore, the plasmon dispersion at low q is 
indeed a fingerprint for the degree of isolation. The strictly 
linear dispersion of the tt plasmon is only visible in isolated 
tubes. 

Although this dispersion might resemble the linear disper- 
sion of 'Dirac electrons' in graphene, the structures seen in 
our VA-SWNT measurements (4-9 eV) are clearly outside the 
energy range of linear band dispersion in graphene. Neverthe- 
less, we will show in the following that graphene is indeed 
the system to be used for the interpretation of the on-axis tt 
plasmon dispersion of VA-SWNT. Our ab initio calculations 
allow us to isolate unambiguously the features of this proto- 
type system. 

To this aim, the loss function — Ime"^((f, oj) of graphene 
was calculated for different momentum transfers q along the 
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FIG. 2: (color online) (a) loss function of the tt plasmon region at 
equidistant g ranging from 0.1 A"^ (top) to 1.0 A"^ (bottom). Right 
stack: observed (filled symbols) vs. calculated (open symbols) tt 
plasmon dispersion for (b) graphite, (c) bundled SWNT [80 vs. dou- 
ble layer graphene and (d) VA-SWNT vs. graphene. The calculations 
are averaged over the two in-plane directions FK and FM. 



FIG. 3: (color online) (a) loss function of graphene at g = 0.41 A"^ 
along FM calculated from the JDOS (green dots), within the bare 
RPA (green) and in RPA including LFE (red, thick). The latter 
changes significantly, when transitions next to the K point (shaded 
area), are excluded (red, dot-dashed), (b) dispersion of the peaks 
in the loss function for different momenta along FM (solid lines & 
squares) and FK (dotted lines & circles) 



in-plane FM and FK directions, for values of (7 = \q\ rang- 
ing from 0.05 to 0.8 A~\ Starting with the bare RPA (with- 
out LFE), the loss function is determined by the independent- 
particle response function xo- as in that case Im oc Im xq- 
The resulting spectra can hence be interpreted as a sum of in- 
dependent transitions, which are directly related to the band- 
structure. Fig. 3a shows a typical spectrum for q = 0.41 
In the low energy (< 10 eV) region, only transitions between 
the TT and 7r*-band contribute to the spectrum, which consists 
of three peaks in FM direction (thin green line) but only two 
peaks for FK (not shown). In Fig. 3b the corresponding dis- 
persions are depicted (green solid and blue dotted lines, re- 
spectively). The first peak arising from transitions within the 
'Dirac cone' at K starts for the lowest q at 0.5 cV and dis- 
perses linearly to 4.0 cV for the largest q. The second peak, 
only visible for FM is a weaker structure around 4 eV which 
shows almost no dispersion. The last peak starting at 4.0 cV 
shows a quadratic dispersion at small q. It can be attributed to 
transitions near the edge of the Brillouin zone close to M. This 
peak is almost undetectable in the joint density of states (see 
dotted curve in Fig. 3a) but is strongly enhanced by matrix ele- 
ments: already at the independent particle level, bandstructure 
alone is not sufficient to describe the spectrum completely. 

When LFE are included in the calculation one determines 
e^^ = 1 + from the full response function X ~ Xo + Xo^x^ 
where the bare Coulomb interaction v reflects the variation 
of the Hartree potential. The inclusion of this term accounts 
for LFE and changes the results drastically (thick red lines): 
whereas induced microscopic components have only little ef- 
fect on the in-plane excitations in bulk graphite | fl6|] . LFE are 
of major importance for the isolated sheets. Most importantly, 
they completely suppress the linearly dispersing low energy 
structure as well as the very weakly dispersing second peak. 



Instead, the peak starting at 4 eV is blue shifted by about 
0.8 cV and becomes the dominant structure in the spectrum. 
Its dispersion is strongly modified: LFE transform the for- 
merly quadratic dispersion into an almost linear one (red line 
in Fig. 3b). One can understand the LFE as a mixing of transi- 
tions that occurs in the inversions when one solves the screen- 
ing equation for x = (Xo ^ ~ v)~^. Therefore, the result- 
ing spectra should consist of mixtures of the formally distinct 
peaks. This can involve a significant energy range. It is there- 
fore most interesting, to analyze whether the linearly dispers- 
ing low energy peak has considerable influence on the spectra 
including LFE. By choosing which transitions we include in 
Xo, we compare the spectra with and without the contribu- 
tions from the the linear region of the tt -bands around the K 
point (i. e. low energy transitions). In the bare RPA loss func- 
tion Xo this region gives rise to the shaded low energy peak in 
Fig. 3a. Despite the very different energy ranges, the final loss 
function after inclusion of LFE is indeed significantly affected 
by the inclusion (red solid line) or exclusion (red dot-dashed 
line) of these transitions: in the latter case the dominant struc- 
ture is strongly reduced (the integrated intensity decreases by 
more than 30%) and is red shifted by about 0.4 cV. There are, 
hence, considerable contributions from low energy transitions 
in the LFE corrected plasmon response. With the mixing of 
transitions of different energies also the different dispersion 
relations mix; therefore, the resulting almost linear dispersion 
is indeed a superposition of the dispersion of the main struc- 
tures in the spectrum, including that resulting from the 'Dirac 
cone'. 

Figure 2d shows the resulting almost linearly dispersing 
TT plasmon of graphene (open diamonds) together with the 
on-axis tt plasmon of VA-SWNT (filled diamonds). The 
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graphene plasmon reproduces qualitatively, and even quantita- 
tively, the experimental findings on individualized VA-S WNT. 
Since our results are completely parameter free, we can con- 
clude that beyond qualitative arguments concerning the tight 
relation of bandstructures and a similar mechanism of the 
LFE, graphene can be studied in order to get insight and quan- 
titative information about VA-SWNT, and vice versa. 

We now move on to elucidate the interacting case of bun- 
dled tubes as shown in Fig. 2c. Experiments on bulk SWNT 
(filled circles Jit]) revealed an asymptotic dispersion relation; 
the TT plasmon is initially shifted to higher energies at small 
q before it approaches the linear dispersion of VA-SWNT at 
large q ^ 0.5 We employ bilayer graphene as an ap- 

propriate model system for representing a typical next neigh- 
bor situation. Indeed, the calculated tt plasmon dispersion of 
isolated double layer graphene with an interlayer distance of 
d ~ 3.3 A (as in graphite) shows the same overall behav- 
ior (open circles) as the measurements on bundled tubes: we 
find a transition from a regime, where the Coulomb interaction 
V oc q~^ yields long range contributions involving neighbor- 
ing layers, to a situation where q is sufficiently large to confine 
main interactions to the same plane. In particular, we studied 
at which interlayer distance d the crossover from interacting to 
non-interacting sheets occurs. For small q = 0.1 A \ a dis- 
tance of 30 A is necessary in order to suppress the influence of 
neighboring sheets on the spectra, while for q = 1.3 A" ^ the 
interlayer distance can be reduced to 7 A. In close analogy, 
the distinct tt plasmon dispersion of bundled SWNT repre- 
sents a smooth dimensional crossover from three dimensional 
bundles to one dimensional separated wires. Hence high q 
measurements are applicable to probe the intrinsic properties 
of individual tubes in bundles. 

Summarizing, we observe distinct tt plasmon dispersions in 
bulk graphite, bundled SWNT and individualized VA-SWNT. 
Owing to the individualization in the VA-SWNT we find a lo- 
calized perpendicular and a strictly linear on-axis tt plasmon 
dispersion. Our ab initio studies uncover drastic changes of 
the spectral RPA response of graphene upon the inclusion of 
local field effects. These LFE account for a linearly dispersing 
TT plasmon in isolated graphene. If a system can be considered 
to be isolated or not depends strongly upon the momentum 
transfer q. In bundled SWNT a transition from an interacting 
to a quasi non-interacting regime for large q occurs and leads 
to an asymptotic dispersion relation. Measurements on VA- 
SWNT assisted by calculations on graphene-based systems 
can hence discern the contributions of the building blocks and 
their interaction, and show that the study of a prototype sys- 
tem of this kind can be used to obtain insight into the collec- 
tive electronic excitations of related materials. 
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